Phonon anomalies due to stripe collective modes in high T c cuprates 
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Phonon anomalies observed in various high T c cuprates by neutron experiments are analyzed 
theoretically in terms of the stripe concept. The phonon self-energy correction is evaluated by 
taking into account the charge collective modes of stripes, giving rise to dispersion gap, or kink and 
shadow phonon modes at twice the wave number of spin stripe. These features coincide precisely 
with observations. The gapped branches of the phonon are found to be in-phase and out-of-phase 
oscillations relative to the charge collective mode. 



The pairing mechanism of high T c superconductors 
has been still highly debated since its discovery over a 
decade ago. The debates are mainly centered on how 
to describe the normal state above its superconducting 
transition T c . Recently, a remarkable series of neutron 
experiments on various cuprates has been done to re- 
veal unequivocally the phonon dispersion anomalies at 
the particular reciprocal point. The phonon anomalies 
differ slightly in their appearance, ranging from a sharp 
dispersion jump to the dispersion kink. This depends 
on the materials; La 2 _ x Sr x Cu0 4 (LSCO)@, 1 1 or 
YBa 2 Cu 3 6+a: (YBCO)§, |, |, § or the focused phonon 
modes such as breathing mode, or bond-stretching mode, 
etc. This ubiquitous phenomenon for which proper ex- 
planation has not been given so far is interesting because 
not only it gives a clue to pairing mechanism in high T c 
cuprates, but more importantly it gives a key to better 
understanding the normal state properties with full of 
mysteries. 

It is known now that in undcrdoping La2- K Sr :E Cu04 
the incommensurate static magnetic peaks at Q are ob- 
served around (tt, 7t) point in two dimensional Cu02 
plane (here the lattice constant a is unity) by neutron 
experiments^. As doping proceeds with x, the order- 
ing vector Q rotates precisely at the insulator-metal 
transition (a; = 0.05) by 45° around (tt, tt) point, cor- 
responding to spin modulation change from the diago- 
nal stripe Q = (27r(i ± S),2tt(^ ± <5)) to the vertical 
stripe Q = (2tt(± ± 5), 0) or (0, 2tt(± ± 5)) where 5 ~ x. 
In higher doping region the incommesurability 8 almost 
saturates at a particular value S ~ i. The same 45° 
rotation of the magnetic superspots is also observed |9| 
in (La, Nd^-zSr^CuO^ corresponding to the metal- 
insulator transition, further lending this superspot ro- 
tation to the universal mechanism to drive a metal- 
insulator transition upon doping. In this material |lC| and 
La2_z(Sr,Ba) 3; Cu04|llJ] the concomitant static charge 
superspots at 2Q are found by neutron scatterings. 

These characteristics coincide nicely with the incom- 
mensurate spin density wave or stripe picture [jl2], [f| 
where a solitonic midgap state accommodates excess 
holes completely (partially) in an insulating (metallic) 
state, giving rise to a microscopic charge and spin seg- 
regation characteristic to topological doping. The par- 



tially occupied midgap state is observed directly by re- 
cent angle- resolved photoemission experiment ]l4|] near 
(tt, 0) point as a parallel Fermi surface sheet which is 
responsible to metallicity of doped CuC>2 plane Q. 

As for YBCO, although experimental verification of 
true static magnetic order comes only from /iSRfTrj] 
which reports a similar phase diagram to LSCO in dop- 
ing vs. static magnetic ordering, there are several com- 
mon characteristics mentioned above |l7||, including neu- 
tron scatterings [[l8|. Thus it is expected that the stripe 
picture may be applicable to YBCO . Numerous other 
experiments such as wipe-out effect for charge order by 
NMR [p0[, d etailed transport measurements J2l), etc. and 
theoriesp2[ , collectively point to support the idea that 
the stripe concept may be applicable to doped Cu0 2 
plane in general. 

Here we are going to demonstrate that the phonon 
anomalies observed in various cuprates are a direct man- 
ifestation of existence of the low-lying collective charge 
excitations associated with stripes. Physically, the tt- 
phase shifting domain walls accommodate excess holes 
and the remaining antiferromagnetic regions are intact. 
The domain wall motion, or the charge collective oscil- 
lation, is strongly coupled with the particular phonon 
mode, such as breathing or bond-stretching modes. Both 
electron and phonon systems are mutually influenced by 
the others to lower the total energy. The spin and charge 
inhomogeneity is a universal phenomenon in doped Cu0 2 
plane, out of which high T c superconductivity is created. 
This process is still beyond our understanding at the mo- 
ment, but this identification, we believe, is important to- 
wards this step. This is particular true for YBCO where 
there is no defining experiment to uncover the charge 
collective modes which has been elusive so far. Since 
the stripe concept, although originally derived from mean 
field approximation may be valid for describing 
a doped Mott insulator, or at least, the stripe is a conve- 
nient "vehicle" to accommodate excess holes, it is worth- 
while to pursue this identification. 

We start with the Hubbard model H = £j 3 - )<r ij ^cj a Cj, a 
+J7Sj7ijjnij to describe a stable stripe and its ex- 
citations. The mean- field ground state jl2], |l5| and 
random phase approximation (RPA)|23| yield a rich 
structure for individual and collective excitation spec- 
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FIG. 1: Dispersion curves of the charge collective mode 
around 2Q. We plot Imx„ n (q,w) for two paths along 2Q 
(a) and perpendicular to 2Q (b). 



tra of spin and charge channels; namely the dynamical 
transverse (longitudinal) susceptibility Xxxizz)^ w) = 
(( i ?i(z)i'Si(z)})q,u! an d the dynamical charge suscep- 
tibility x„„(q, uj) = ((n; n)) q>(i ,. The RPA equation 
for xs^ = ((5 , T; S 'i>> is written as xs T Sj (ri, r 3 , w) 
= Xo J, ( r i> r 3><*>) + C/ Er 2 Xo i ( r i> r 2j w ) XS T Si (r 2 ,r 3 ,u;) 
with non-interacting susceptibility xj^- After Fourier 
transformation to k-space, this is reduced to a matrix 
equation of xj an d XSfSi- By solving it, we obtain 
Xxx(q + hQ, q + hQ, w). There, q + hQ and q + hQ 
are, respectively, outgoing and incoming wave vectors. 
They can change by (h — h)Q through the Umklapp 
process, since the underlying electronic state is spatially 
modulated with the periodicity Q. In a similar man- 
ner, we calculate ((nt;n|)) and ((nt;ni)), and obtain 

x^ z (q + ^iQ,q + hQ,u) and Xnn(q + 'iQ>q+^Q,w). 

The spin modes Xxx(zz){<\, u) = Xxx(zz){<\,^u) give 
rise to collective excitation at the ordering vector Q 
and its odd-harmonic points as a Goldstone mode. The 
charge mode Xnn(q, w) = Xrm(q, q, exhibits the collec- 
tive excitations at q = 2Q, 4Q • • •. There, Xxx gives spin 
wave, when we assume that the magnetization points to 
the z axis. The translational phason modes appear in 
Xrm iXzz), which correspond to the motion of the charge 
(spin) stripe p^|. Among various individual and collective 
motions associated with stripes, the phason mode in the 
charge oscillation is directly coupled to the phonons and 
most relevant to the phonon renormalization. In Fig.l, 
ImXnn(q, w) is shown around q = 2Q = (0, ?) along the 
q y - and the (^-directions. This particular charge mode 
with the anisotropic velocities is situated at 2Q. The 
motion parallel (perpendicular) to 2Q vector corresponds 
to the compression (meandering) motion of stripes with 
high (low) velocity. In the following we illustrate typ- 
ical results for the metallic vertical stripe |H| which is 
stable with Q = (ir,^-) for the electron filling per site 
n = 0.82 and the parameters are fixed as y = 4.0 and 

j = —0.2 (t(t') is the (next) nearest neighbor hopping 
integral). The main features of the results do not much 
depend on the choice of these values. We note that t is 
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FIG. 2: Renormalized phonon spectral functions for \g\ 2 = 
0.001. (a) is logarithmic plot to emphasize the small inten- 
sity of the shadow bands, (b) is contour plot shown for high 
intensity. The simple sinusoidal form is assumed for the un- 
perturbed phonon (A = -0.01, B = 0.08, k = 0)Bj). 
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FIG. 3: Schematic figure for relevant modes. 1: Charge col- 
lective mode. 2: Unperturbed phonon. 3,4, and 5: Phonons 
shifted by 2Q, 4Q and 6Q via Umklapp processes. 



an order of 0.5eV and the phonon energy uj ~ 80meV 
for the breathing mode or bond-stretching mode which 
is coupled strongly to the electronic state. We focus on 
an energy range of typically j ~ 0.1 in the following. All 
the energy is scaled by t. 

In order to study the effects of the stripes on the 
phonon Green's function D(q+ZiQ, q+^Qi * w n) within 
RPA, we consider the following standard electron- phonon 
Hamiltonian: H' = Sk, q ,<rffCk +q!0 .c k!(T (6 q + Z^ q ) with g 
being the coupling constant. Since the electron bubble 
in the phonon self-energy correction is nothing but the 
dynamical charge susceptibility mentioned, we immedi- 
ately obtain an equation in a matrix form for the thermal 
Green's functions £ m '{<W - |ff| 2 A)(q+^Q, iw n )x™(q+ 
IQ, q + m'Q, zcj„)}L>(q + rn'Q, q + rnQ, iu n ) = Do(q + 
lQ,iui n )5i m , where the unperturbed Green's function 
£>o(q, iu> n ) = — 2w q / (w^+w q ), and w q is the unperturbed 
phonon dispersion. Diagonalization of this matrix equa- 
tion yields the renormalized phonon Green's function. 

We show the renormalized phonon spectral function 
— ilm£)(q, q, iu> n — > uj + iff) in Fig. 2, where we assume 



a simple sinusoidal form|24| for w q as a typical phonon. 
It is seen from Fig. 2 that the dispersion exhibits several 
characteristic modifications: (A) The discontinuities of 
the dispersion occur at various points, notably near 2Q 
where the dispersion cone of the collective charge mode 
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FIG. 4: Schematic figure for two in-phase (lower energy) and 
out of phase (higher energy) oscillation patterns correspond- 
ing to the two branches above and below the gap at 2Q. The 
vertical bars show the charge domain walls of the stripes. 



shown in Fig.l and the phonon mode intersect. (B) The 
appearance of the downward hybridized new phonon dis- 
persions. (C) The flat dispersion appears above the gap 
at 2Q. (D) The shadow band which is a replica of the 
original mode is seen. 

These features come from the two physical reasons 
schematically illustrated, in Fig. 3: The hybridization of 
the charge collective mode centered at 2Q denoted as 1 
with the original phonon 2 in Fig. 3 and the band-folding 
due to the periodic structure with the charge ordering 
vector 2Q, which amounts to shifting the original dis- 
persion 2, resulting in phonon dispersions; curve 3 ~5 in 
Fig. 3. The hybridization of these dispersions gives rise 
to the reconnected phonon modes in Fig. 2. 

These main four features (A-D) of the dispersion mod- 
ification become more evident when the electron-phonon 
coupling \g\ 2 increases, or the unperturbed dispersion is 
situated at lower energy side. The gap size at the discon- 
tinuity increases as the original phonon band is widen. It 
is noteworthy that the discontinuity of the dispersion, or 
a kink feature and associated new lower mode branching- 
out from the disconnected lower phonon are reminiscent 
of the observation on 06. 35 and Og.60 in YBCO by Mook 
and DoganQ (see their Figs. 1(d) and (e)). Note also that 
since the discontinuity point occurs at the intersection of 
the charge collective mode at 2Q and phonons, the ob- 
served dispersion anomaly should change with doping; As 
doping increases the 2Q vector increases. This is indeed 
observed by neutron experiment |6) and also consistent 
with the doping dependence of the fundamental vector 
Q of the magnetic peaks These constitute an inter- 
nal consistent picture for the stripe origin. 

The physical implication of the phonon gap formation 
is understood by analyzing the phonon oscillation pat- 
terns of the upper and lower phonon branches. Accord- 
ing to linear response theory, the oscillation patterns are 
given by (f>(r,t) = £^£>(q + ZQ, q, u)e il ^ T h tl e i ^ T - i " t 
where </>(r, t) is the lattice displacement and /i q is the ex- 
ternal field with the wave number q and the frequency u. 
At the wave number 2Q = (0, ^), the oscillations corre- 
sponding to the upper (w = 0.08) and lower (u = 0.075) 
branches shown in Fig. 2 are schematically displayed in 
Fig. 4. The four site period oscillations are either in phase 
or out of phase with the charge oscillation in the electron 
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FIG. 5: Renormalized phonon spectral functions for two 
perpendicular paths. \g\ 2 = 0.001. (a) and (c) (0, q y ) for 
< q y < 7T. (b) and (d) (q x , § ) for \q x \ < f. The unper- 
turbed phonon with a steep slope (A = 0.01, B = 0.08, k = 
0.9999) |24| is assumed, (a) and (b) are logarithmic plots, (c) 
and (d) are contour plots. 



system which is the phason motion of the stripes. 

We show the renormalized dispersions parallel to 2Q 
in Figs. 5(a) and (c) and perpendicular to 2Q centered 
at 2Q in Figs. 5(b) and (d). Here we assume the unper- 
turbed phonon with a steep slope at (0, -|), simulating 
the breathing mode in LSCO with the energy ^80meV. 
Noticeable features of the modification are: (A) Rather 
large dispersion discontinuity observed at 2Q is seen in 
Figs. 5(a) and (c). (B) It is noted that the shadow mode 
which is a replica of the original band can be seen in 
Figs. 5(a) and (b). (C) The anomaly continues to ex- 
tend to the perpendicular direction seen from Figs. 5(b) 
and (d). These features (A-C) coincide with the neutron 
experiment^] on LSCO with x = 0.15 where at (f,0) in 
the (q x ,0) scan they observe a sharp drop of the bond- 
stretching mode (see their Fig. 2) and this anomalous be- 
havior continues to the perpendicular direction (see their 
Fig. 4). As they emphasize, the phonon anomaly occurs 
not because of the simple band-folding effect due to new 
charge periodicity with 4 lattice period which would re- 
sult in a gap at the new Brillouin zone boundary at (f, 0), 
but the observed gap is at (f ,0). 

The anomalous reciprocal regions for the two perpen- 
dicular directions have different widths as seen from Fig. 5 
because the underlying charge cone is anisotropic as men- 
tioned, yielding different widths (compare the openings 
of the cones at the center in Figs. 1(a) and (b) correspond- 
ing to the cones in Figs. 5(a) and (b)). Thus it is possible 
to use it as a kind of spectroscopy, namely, to investigate 
this anisotropy of the charge collective mode through the 
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analysis of the phonon anomaly. In fact according to the 
above experiment Q the width of the anomalous region 
in the perpendicular direction is wider than that in the 
parallel direction. This accords with our result in Fig. 5. 
We can attribute it as arising from the lower velocity of 
the meandering mode than the compression modep3|. 

McQueeney, et al.fn report that there is no phonon 
anomaly below x < 0.05 in LSCO where the static di- 
agonal stripe is stabilized. There are two possibilities: 
(1) Since the velocity of the collective charge excitation 
is rather low for the insulating diagonal stripe |23|, the 
phonon mode investigated by them may not be coupled 
with this excitation, failing to drive the anomaly. (2) 
As the second possibility, the collective charge excitation 
could be gapped |^| . There is no low-lying excitation 
coupled to phonons. In this connection, we point out 
that in Cr, where the diagonal stripe in the present termi- 
nology is realized, neutron experiment to see the phonon 
anomaly at 2Q corresponding to the so-called strain wave 
or charge order is desirable. 

Considering the electron-phonon coupling, we study 
the phonon anomalies in the stripe state. The phonon 
anomalies observed mainly by neutron scattering exper- 
iments on LSCO and YBCO with various doping levels 
are accounted for as arising from two physical effects: 
Hybridization with the charge collective mode at 2Q and 
band-folding of a phonon band. These give rise to re- 
connection and gap formation in bands. The underly- 
ing physics behind the phonon anomalies at the partic- 
ular reciprocal point, namely 2Q, is the existence of the 
collective charge excitations associated with stripes; the 
translational phason and meandering gapless modes with 
different velocities situated at every even-multiples of the 
fundamental vector Q which varies systematically with 
hole-doping. Thus this also leads to systematic changes 
of the 2Q vector, a fact that the experiments have found 
(see Fig. 4 in Ref. |j|). In order to confirm our identifica- 
tion, further neutron experiments are desirable, in partic- 
ular to find the yet-uncovered lower modes branching-out 
downwardly from the phonon discontinuity at 2Q. Since 
in principle the phonon anomalies occur for all phonon 
modes and the observability depends on the electron- 
phonon coupling constant, lower phonon modes are more 
favorable for this purpose. 

It might be interesting to point out that similar phonon 
mode anomalies are also found in (La,Sr)2Ni04(25) where 
the diagonal stripes are stabilized because as mentioned 
this system is an insulator up to higher dopings, thus 
45° rotation of the superspots occurs. In yet another 
non-cuprate superconductor Bao.6Ko.4t3i03, which is re- 
lated to charge density waves instead of spin density wave 
here the phonon anomaly of the bond-stretching mode 



appears at a certain reciprocal point j2(|. We believe that 
these anomalies could be analyzed in the same way as in 
here, which belongs to future problem. 

We thank J. Tranquada, R. McQueeney H. Mook and 
N. Wakabayashi for their useful discussions and informa- 
tion. 
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